Crystal-symmetry preserving Wannier states for fractional chern insulators 
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Recently, many numerical evidences of fractional Chern insulator, i.e. the fractional quantum 
Hall states on lattices, are proposed when a Chern band is partially filled. Some trial wavefunctions 
of fractional Chern insulators can be obtained by mapping the fractional quantum Hall wavefunc- 
tions defined in the continuum onto the lattice through the Wannier state representation (Phys. 
Rev. Lett. 107, 126803 (2011)) in which the single particle Landau orbits in the Landau levels 
are identified with the one dimensional Wannier states of the Chern bands with Chern number 
C = 1. However, this mapping generically breaks the lattice point group symmetry. In this paper, 
we discuss a general approach of modifying the mapping to accommodate the lattice rotational 
symmetry. The wavefunctions constructed through this modified mapping should serve as better 
trial wavefunctions to compare with the numerics and also as the basis for construction of lattice 
symmetry preserving pseudo-potential formalism for fractional Chern insulators. The focus of this 
paper shall be mainly on the C4 rotational symmetry of square lattices. Similar analysis can be 
straightforwardly generalized to triangular or hexagonal lattices with Cs symmetry. We also gener- 
alize the discussion to the lattice symmetry of fractional Chern insulators with high Chern number 
bands. 
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I. INTRODUCTION 

Integer and fractional quantum Hall states differ from 
the conventional phase of matter characterized via sym- 
metry breaking by their distinct topological proper- 
ties. Under strong magnetic field, the 2D electron 
gas forms highly degenerate Landau Levels (LL) which 
carry non-trivial topological indices known as the TKNN 
invariants^ or the Chern number. When a Landau level 
is fully occupied, the Chern number reveals itself as the 
quantized Hall conductance in integer quantum Hall ef- 
fect (IQH effect). When it comes to a partial filling, the 
strong correlation between electrons arise from the huge 
degeneracy of LL (or the quenched kinetic energy), to- 
gether with the its non-trivial topology, leads to the frac- 
tional quantum Hall states (FQH states) that support 
fractionalized excitations and non-trivial braiding statis- 
tics. In a work in 1988^, F. D. M. Haldane proposed the 
first realization of a band insulator model with non-trivial 
Chern bands (i.e. bands with non-zero Chern number). 
A fully occupied Chern band also shows a non-zero Hall 
conductance, which is referred to as quantum anomalous 
Hall (QAH) insulator, or Chern insulator. This resem- 
blance between the Chern bands and that LLs is rooted 
in their equivalence in the sense of topology. Now it is 
conceivable that, if the Chern band is flat or almost flat, 
fractional Chern insulators (FCI), which are FQH states 
without magnetic field, will emerge. Following this line 
of thought, lots of effort has been devoted into the search 
for FCI states3-i9. 

For the FQH states, wavefunctions, the validity of 
which can be checked most straightforwardly through 
numerics, have been a major tool to study the underly- 
ing physics. However, the understanding of FQH states 
through wavefunctions such as the Laughlin wavefunc- 
tion is not directly applicable to FCI states, since the 
former can be written as holomorphic functions and the 
latter is defined on a lattice. As proposed in Ref. [H, this 



difficulty can be overcome through the mapping between 
the LLs and Chern bands that exploit their similarity 
in the Wannier state representation (WSR). In a Chern 
band, the single particle Wannier states that are local- 
ized on one direction and carry well-defined wave vec- 
tors on the other resemble the Landau orbits in the LLs. 
By directly identifying them in the Fock basis, all FQH 
wavefunctions written in the continuum can be mapped 
into FCI wavefunctions. In the WSR, there is a gauge 
choice in the definition of the Wannier states in the Chern 
bands. Different gauge choices, though capturing the 
same topological properties, lead to microscopically dif- 
ferent FCI states through the mapping. Attempts have 
been made to gauge-fix the Wannier states for the com- 
parison with numerics22ri2^. However, these gauge fixing 
schemes generically render the lattice rotation symme- 
try broken in the mapping to the LLs unless the Berry 
curvature in the Brillouin zone is homogeneous^. An ex- 
ception is the Coulomb gauge in Ref. [20 which preserves 
the four-fold lattice rotation symmetry, though it is not 
constructed for this purpose. This gauge choice will be 
further discussed later. 

In this paper, we discuss the general procedure of con- 
structing Wannier states that preserves crystal symme- 
try. Rather than using WSR, we equivalently identify the 
Bloch states of the LLs and Chern bands, similar to Ref. 
[20I . We show that there is a series of gauge choices which 
leads to Bloch states in Chern bands that transform in 
the same fashion as that in the LLs. With these gauge 
choices, coherent states for the Chern bandsi^ preserv- 
ing the lattice rotation symmetry can be constructed. 
This implies that the FCI states and their correspond- 
ing pseudo-potential Hamiltonians constructed through 
the mapping between the Bloch states in the LLs and 
the Chern bands are invariant under the lattice rotation. 
In the following, we will first introduce the Bloch state 
representation of the LLs and discuss its rotation proper- 
ties. Then we will use the lattice Dirac model-- and the 
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Checkerboard model- as examples to demonstrate the 
construction of the rotation symmetry preserving gauge 
for the Chern bands with C = 1. In the end, we will 
discuss the case with higher Chern number. 



II. GENERAL CRITERIA 

In the original proposal^^, a LL is mapped to a Chern 
band by identifying the Landau orbits with the Wannier 
stats. Equivalently, we can think of this mapping as an 
identification between Bloch states of the LL (that will 
be defined later) and those of the C = 1 Chern bandsi^ 
Although not explicit in this Bloch state representation, 
the LL does have a continuous SO{2) rotation symmetry 
about the z axis which includes a lattice rotation symme- 
try as a subgroup. Thus, lattice rotation symmetric FCI 
states can be obtained from 5*0(2) rotation symmetric 
FQH states in the continuum written in the Fock basis 
by mapping the LL to the Chern bands in the Bloch state 
representation and requiring that the Bloch states in the 
LL and the Chern bands share the same transformation 
properties under the lattice rotation. As the first step, 
we should define the Bloch states of LL and study their 
behavior under lattice rotation. In the rest of the pa- 
per, we will mainly focus on the C4 rotation symmetry 
of square lattices. 

In the Landau gauge, the Landau orbits can be written 

as 
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where Ly is the size of the system in y direction, li, = 
yjhj (eB) is the cyclotron length and Py = 2-Kn/Ly, n € 
Z is the wavevector along y direction. As we can see, 
the Landau orbits are plane waves along the y direction 
but localized along the x direction. Now, we introduce 
a square lattice to the system with lattice constant a 
satisfying = 27rZ^. Since each unit cell contains one 
fiux quantum, the translations along the two base vectors 
commute with each other, which allows us to define Bloch 
states as an alternative basis of the LL. Consider an x 
N system (in terms of lattice constant a). The Bloch 
wavefunction can be written as 

<(x, y) = ^ e''^--"''V'fc„+2„v<^(^' yl (2) 

where k = (fc^,, ky) is the Bloch momentum with kx,ky S 
[—■^, f ] ^^'^ ^^.y ~ '^T^iT'x^y/N, rix^y G Z. Note that 



(3) 



where g^^j, are the two base vectors of the reciprocal lat- 
tice. The Bloch wavefunction can be expressed in the 
standard form 



<(x,y)=e^'=-^+^'=«^Wk(x,2/), 



(4) 



with 



Mk = 



E ikx {7ia—x)-\'i- 



(5) 



One can show straightforwardly that the Berry connec- 
tion A in the Brillouin zone (BZ) is given by 



A = -i(uk|Vk|uk! 



= i-lf.k 



0), 



(6) 



from which we see that the Chern number of the LL 
is C = 1 by integrating the Berry curvature across the 
whole BZ. From the definition of Uk{x,y) in Eq. [SJ one 
can derive that 



Uk{y,-x)^e i^-'^y'-le Uflk (a;, y). 



(7) 



where i?k — {—ky, kx) is the Bloch momentum k rotated 
by 90°. Eq. [7| describes how the Bloch wavefunctions 
transform under the C4 rotation in real space. On the 
R.H.S. of Eq. [71 uiik(x,y) is what would be normally 
expected for Bloch states under C4 rotation. The fac- 
tor e~^^yl^b generates a real space gauge transformation 
which follows a spatial C4 rotation and restores the ap- 
parently broken rotation symmetry by the Landau gauge. 
The factor eT'^^'^^y^^ results from a specific gauge choice 
of the Bloch states in the BZ of the LL. 

To extract the gauge condition of the Bloch states that 
is crucial for C4 rotation symmetry of the LLs, we will 
introduce the coherent states i/izj, of the LL which form 
an overcomplete basis: 



'^zo{x,y) 
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where zq = xq + iyo is the center of mass of the coher- 
ent states. These states are localized in both directions. 
And they are invariant under the 50(2) spatial rotation 
around its center up to a real space gauge transformation. 
We can re-express them in terms of the Bloch states: 
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(9) 



This expression will serve as the general definition for 
coherent states on Chern bands when they are mapped 
to the LLs. We notice that it is the gauge condition 
Eq. [3] that guarantees the localization of -020 • One also 
can verify that, under the 04 rotation, the coefficients 
in the second line conspire with the factor e"**^"'*^"'' to 
guarantee the invariance of ipzQ under 04. 

Thus, if we want to establish an identification be- 
tween the Bloch states of the LL and those of the Chern 
bands without breaking the lattice rotation symmetry, 
the Bloch states of the Chern bands should satisfy the 
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same gauge condition: (1) Bloch states has the same pe- 
riodicity as that in Eq. [31 (2) Under C4 rotation, the 
Bloch wavefunction Uk satisfies: 



(10) 
(11) 



where Rc^ denotes the symmetry transformation inside 
the unit cell, which is (x^y) — >■ {—y,x) followed by a 
real space gauge transformation in the LL case. As a 
convention, we can always require that Rq^ — 1 and 

Uk=o is the eigenstate of Rc^ with eigenvalue 1. 

III. EXAMPLES 

Lattice Dirac model - Now we use the lattice Dirac 
model as the first example to demonstrate how to pre- 
serve the C4 symmetry by a gauge choice. The lat- 
tice Dirac model is given by the Hamiltonian H = 
X]k Ck'H(k)ck, where Ck is a two-component fermion op- 
erator and 



H(k)=a-dk (12) 



with 



dk = (sin kx, sin ky, M + cos k^ + cos ky). (13) 

We have taken the lattice constant to be a = 1. M is a 
free parameter which is taken to be < M < 2 in order 
that the lower band has C = 1. This model is symmetric 
under a simultaneous C4 rotation in both the real space 
and spin space: 

'H( ky, kx) = Rc^Hikx ■, ^y)^Ci ' 



(14) 



where Rc^ — e 







[Ij ) is the spin rotation ma- 



trix times a U{1) charge rotation which is choosen, for 
later convenience, to satisfy i?^^ = 1. 

To gauge fix this Chern band, we can take the following 
two-step strategy. First, we write down the Bloch states 
in a gauge that respects the full periodicity of BZ and 
the C4 rotation symmetry. Usually, this gauge will have 
singularities at the high symmetry points of the BZ. For 
the second step, we modify the gauge to satisfy the gauge 
condition Eq. [TUland Eq. [TTJ Since Eq. [11] only connects 
the Bloch states with momentums that are related to 
each other by C4 rotation, we can pick any gauge in the 
"principle region" (PR) which is defined as a quarter of 
the BZ, kx^y e [0, 7r/a] as shown in Fig. [T](a), in this case. 
Then we use Eq. [Tl]to generate the gauge on other parts 
of the BZ. The requirement of self-consistency and the 
consistency with gauge condition Eq. [TUj imposes only a 
boundary condition to the the gauge choice in the PR. 

Now, for the lattice Dirac model, we first write down 
the Bloch states in the lower band: 



kx , k^. 



(15) 



where cos^k = d^/|dk| and Lpk = arg((i^ -f zdk)- This 
gauge satisfies 

\kxTky) — \kx 27r, ky) — \kx , ky ^ 2tt) , 

\ - ky,kx) = Rc^\kx,ky). (16) 

It is easily shown that \kx,ky) is singular only at k = 
(tt, tt). For step two, we denote the Bloch states in the 
Dirac model which maps to those of the LL with the C4 



symmetry preserved as \kx,ky)) = e 



i4>(kj: ,ky) 



\kx,ky). The 



gauge conditions Eq. [TUj and Eq. [TT] then requires that 



^i4>{k^,ky) _ ^i4>{k^+2TT,ky) ^-^j-^ 

Notice that the "twisted" periodicity gauge condition 
along the ky direction in Eq. [TUjcan be generated by these 
two equations. Therefore, the phase factor e^'^'^'^"'^^^ in 
the BZ is completely defined by its value in the PR with 
the following consistency condition: 



(^(fc,7r) 



.^(fc,0) =0(O,fc, 
k _ 
2 ~ 



vr, k) = (/>(7r, k), 



(18) 



where k € [0, tt]. Notice that the second equation leads 
to a singularity of 0(k) on the boundary of the PR, since 



lim (/)(fc, tt) + 7r/2 = lim (/>(7r, k). 

k^TT k^TT 



(19) 



This singularity correctly cancels that of fcj^) at kjvi = 
(tt, tt). All function 0(k) in PR that satisfies the gauge 
condition Eq. [TSj and that are smooth only except at km 
can be consistently extended to the whole BZ by Eq. [TT] 
As an example, one solution is given by: 



< 1 



TT ; k' TT ^ ^ (20) 

0, elsewrhere in PR 



where 9'^ = aTg{kx + iky — (1 -I- «)7r) — tt. Now we have 
obtain \kx,ky)) which corresponds to Uk in the LL. With 
this identification, coherent states for the Chern band in 
the lattice Dirac model can be defined through Eq. [9l 
and the C4 rotation properties will be inherited from the 
LL. 

At this stage, it seems that the cancelation of singular- 
ities between 0(k) and \kx, ky) is quite accidental. How- 
ever, it is shown in Ref. [25] that the symmetry enforces 
constraints on the Chern number. And we will show 
that the Chern number dictates the singularity and leads 
to this cancelation. In the first step, we have chosen a 
gauge that satisfies Eq. [THj and usually leads to singu- 
larities at the high symmetry points: the F point with 
kr = (0,0), the X point with kx = (ti'iO), the Y point 
with ky = (0,7r) and the AI point with k^/ = (7r,7r). 
The singularities at these points are U{1) monodromy 
which can be determined by the eigenvalues of i?c4 
these points. At the F point, by definition. 



lim Ikr 

e-)-0 



egy) = RcJ'kr + egx)- 



(21) 
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where 

, kx ky . kx . ky cos kx cos ky 
h, = (cos - COS - sm - sm ) 

(25) 

and ek — ^ ^ cos kx cos ky that does not affect the Bloch 

\/2+2 " 

states. The BZ is defined 1 with the 

gauge condition 



FIG. 1. (a) The BZ of square iattice with kx,y £ [— 7r,7r] is 
indicated by the dash fine. The biue area is the PR with its 4 
corners being high symmetry points: F, X, Y and M. (b) The 
red fine is a contour that foiiows the boundary of the PR. ab, 
cd, ef and gh are quarter circies with inhnitesimai radii that 
centered around the 4 corners of the PR. 



<J^T-L{kx,ky)a^ = H(fc^ + 27r, ky) = T-L{kx, ky + 2ii). 

(26) 

The C4 symmetry and gauge property of this model is 
given by 



a'x'H,{ kyjkx^cFx — 'H-{kx,ky). 



(27) 



Meanwhile, |kr) must be an eigenstate of the with 
its eigenvalue denoted as ^r- Thus, we have 



|kr + egy) -Cr|kr + egy 



0. 



(22) 



Thus, the U{1) phase singularity along the infinitesimal 
quarter circle ab is ^p^, and so is the Berry connection 

integral along the path ab in this gauge. Parallel analysis 
can also be done for M point relating the inverse of the 
eigenvalue of Rc^ at the M point denoted as £,m with 
the U{1) phase singularity and also the Berry connection 
integral along the path ef. For X and Y points, they 
are only invariant under Rq^- We denote the eigenvalue 

of Rq^ at X as rjx which equals that at the Y point. 
Through similar reasoning, the 773^^ equals the total U{1) 
phase singularity and the Berry connection integral along 
the path cd and gh. Now, we consider the total Berry 
connection integral along the red contour in Fig. [T] (b) 
which should be, when raised to the exponent, e*^^, 
namely a quarter of the total Berry curvature (on the 
exponent). Due to the gauge condition Eq. [111 the Berry 
connection integral along ha cancels that along be. Same 
cancelation happens for de and fg. Thus, the only con- 
tribution to the integral comes from the four infinitesimal 
quarter circles: 



(23) 



This result is also obtained in Ref. l25|. As explained 
above, ^CjT/'yx^ ^^^^ characterized the total U (1) phase 
singularity along the four quarter circles at the 4 corners 
of the PR. For C = 1, comparing this total singularity 
of the gauge \kx, ky) with singularity of the phase factor 
g«0(k) required by Eq. [THl we see that they cancel each 
other. 

Checkerboard model - As a second example system 
where C4 symmetry acts differently, we consider the 
Checkerboard model that produces a flat Chern band^: 



With the reference of the gauge condition Eq. [TU] and 
Eq. [11] we need to find the Bloch states in a gauge such 
that: 



■2TT,ky)) =cr^\kx,ky)) 

- ky,kx)) — —cFx^^ 2 \kxTky)). 



(28) 



Following the two-step analysis described in the lattice 
Dirac model case, we obtain the general solution: 

|fc.,/c,))=e'^('')(j;|j:j, (29) 

where cos 0k = ^k/l'^kl, '■P\l = arg(^k + *^k) ^^"^ e^'^^^'^ is 
a J7(l) phase that satisfies the following condition only 
on the boundary of the PR: 



0(fc,O)=(/.(O,fc) 

'^(I'.fc) = ^ f + 0(fc,7r), 



(30) 



H(k) =ek/ + /i£tT^ 



(24) 



where fc G [0,7r]. As expected, the singularity of the 
(/)(k) around the boundary of PR can be set to cancel 
the singularity of the two-component spinor in Eq. 1291 
Now C4 symmetric coherent states can be constructed 
via replacing ip^^ in Eq. [Hi with the Bloch states in Eq. 
[29l The wavefunction of a coherent state that centers at 
the origin is depicted in Fig. [2](b). 

In the end of this section, let's discuss why the 
Coulomb gauge in Ref. [l^ preserves the C4 symmetry. 
In the Coulomb gauge, the gauge potential Ak is the 
sum of two parts A = Aq -t- SA with Aq identical to the 
gauge potential in LL, and the correction SA with zero 
Chern number. Once gauge fixed by the Coulomb gauge 
condition, SA is single-valued in the BZ and transforms 
trivially under C4 rotation. Therefore the C4 rotation 
property of the Coulomb gauge is completely determined 
by Aq and thus is identical to the LL. Therefore the 
Coulomb gauge preserves the C4 rotation symmetry in 
the mapping to the LL. 
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FIG. 2. (a) On the checkerboard lattice, the red and blue sites 
are the upper and lower components of the Bloch state. The 
coordinates are given in the unit of lattice constant, (b) The 
wavefunction of a coherent state that centers at the origin is 
symmetric under C4 rotation. 



IV. GENERALIZATION TO HIGHER CHERN 
NUMBER SYSTEMS 

So far we have demonstrated that, for C — 1 Chern 
bands, we can always find the appropriate gauge such 
that the C4 symmetry is preserved when the Chern band 
is mapped to the LL. Now we generalize the discussion 
to the case of Chern bands with Chern number C > 1. 
It has been shown in Ref. H) that a Chern band with 
Chern number C can be mapped into C layers of Chern 
number 1 bands by WSR (given the condition that the 
system size is commensurate with C along the direction 
in which the Wannier states are localized). An alter- 
native approach has been developed in Ref. [l^. Since 
this mapping generically breaks the rotation symmetry, 
the existence of a rotationally symmetric scheme requires 
further inspection. In the following we will focus on the 
C4 symmetry of FCI's in a C = 2 Chern band. 

As is discussed in Ref. [H, in general, one-dimensional 
Wannier states of the C = 2 band can be defined by 
one-dimensional fourier transform of the Bloch states |k) 
along a reciprocal vector direction nx+my. For each pair 
of mutually prime integers (n, m), a set of Wannier states 
are defined, which defines a one-to-one mapping between 
bilayer FQH states and the FCI in this C — 2 band. 
If n is odd, the resulting FCI states is sensitive to the 
translation along x direction, such that an x-dislocation 
with Burgers vector h — x will generically become a non- 
Abelian defect. Similarly, if m is odd the y-dislocation 
with Burgers vector h = y becomes a non-Abelian de- 
fect. Therefore the topological properties of the FCI state 
necessarily break C4 symmetry unless n and m are both 
odd. In the following we will study the C4 symmetry of 
the simplest case (n,m) = (1,1) and discuss the more 
general cases later. 

The Wannier states for (n,m) = (1,1) are defined as 



ky) 



(31) 



fc+e[0,27r) 



with L = Lx = Ly the linear size of the system along 
both directions, and Ik^, ky) the Bloch state of the C = 2 



band. k+ denotes the momentum along the diagonal di- 
rection and ky denotes the momentum along y direction 
at fc+ — 0. The Bloch state can have different gauge 
choices and the goal of the following discussion is to find 
the conditions the C4 symmetry impose to the gauge 
choice. For C — 2 bands, the Wannier states at even and 
odd sites form two families which are not connected by 
adiabatically changing ky. Each family is topologically 
equivalent to a C = 1 band, or a Landau level. Since the 
FCI state constructed by WSR are bilayer FQH states 
with each family mapped to a layer of Landau level, the 
C4 symmetry will be preserved if it is preserved for each 
layer. The Bloch states of each family can be obtained 
by an inverse Fourier transformation to each family of 
Wannier states in Eq. [3TJ 



\k+,k+ +ky)^ = E \^2n-l,ky) ' 



|fc+, fc+ + ky)^ 



-2ink-L 



(32) 



Here fc+ e [0,7r) and ky = [0,27r) defines the reduced 
Brillouin zone. Combining Eq. [31] and [32l we obtain the 
following simple relation between the Bloch states of the 
C = 2 states and that of the effective C = 1 bands: 



|fc+, fc+ + ky)-^ = — r=- {\k+,k+ + ky) - |fc+ + 7r,k+ + ky + tt)) 
v2 

\k+,k+ + ky)^ = —= (|fc+, fc+ + ky) + \k+ + TT,k+ + ky + n)) 
v2 

or equivalently, in the conventional kx,ky basis, 



|k), = -^(|k)-|k+(7r,7r)))e-^-^ 



Ik),- 



1 

1 
71 



(|k) + |k+(7r,^))). 



(33) 



To preserve C4 symmetry in the two C = 1 bands, the 
Bloch states |k)j^ ^ need to satisfy the C4 invariance con- 
ditions (Uni) and (dH). Eq. ES] together with Eq. [TUland 
[TT] determines the general C4 symmetry condition for the 
C = 2 FCI states constructed in the (1, 1) WSR. 

For more generic (n, m) with n, m odd, the discussion 
is exactly parallel and the reduced bands |k)^ 2 de- 



[3| with the momentum Q — (tTjTt) 
{mTjimr), and the factor e"**^^ re- 
with s,t two integers satisfying 



det 



= 1. For n,m odd, Q' and Q are equal 



termined by Eq. | 
replaced by Q' = 
placed by e"*'^*''^" 
n m 
s t 

modulo a reciprocal lattice vector. Therefore the effec- 
tive C — 1 bands for all (n, m) odd are equivalent up to 
a phase factor. 

In particular, the C4 invariance condition imposes 
some constraints on the rotation eigenvalues of the C = 2 
band at the high symmetry points. The original and re- 
duced BZ's are shown in Fig. [3] In the same way as the 
C = 1 band discussed in the last section, we denote the 
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k> 



M„ 



FIG. 3. A Chern number C = 2 band is mapped into two 
subbands with C = 1 by folding the BZ, namely taking a 
superposition of the states |k) with |k + ki/o). 



i?C4 eigenvalue of the original C = 2 band at high sym- 
metry points Tq,Mq as ^ToAmq, and the R^^ eigenvalue 
at Xq point as rjxo ■ At the F points of the reduced BZ, 
the C = 1 Bloch states |k) 2 ^^^^ superpositions of | (0, 0)) 
and Ktt, tt)) in the original BZ. Therefore the C = 1 bands 
at F point can only be C4 invariant if — Cmq ■ The 
i?C4 eigenvalue of the two C = 1 bands at the F point 
satisfies ^ri,2 = fro = Cmq- At the M points of the sub- 
bands, the Bloch states are the superpositions of |kxo)o 
and |kY-(,)o which are both eigenstates of Rq^ with eigen- 
value rjxo- Hence, the eigenvalue of Rc^ at Mi ^2 should 
satisfy ^ = rjxo ■ At the X points of the subbands, the 
Bloch states are the superposition of |k = {—^, — f )) and 
|k = (^, ^)). Since Rf^^ = 1 by definition, the eigenvalue 
of R(j^ at Xi^2 should be rjx^ ^ = 1. By construction 
the two bands |k).^ 2 each have Chern number C = 1, so 
that we obtain based on Eq. [221 that ^r^^^i'^Vxl = ^1 
and Cr/^M^^Xj^ ~ * with j = 1,2. To summarize, the 
C4 symmetry requires the following conditions to the ro- 
tation eigenvalues of the original band and the effective 
C — 1 bands: 



^Fo^mIvxI - 1' 

f2 _ /:2 _ 2 _ 2 

— ?Af2 ~ '^^o' '7-^1 ~ ^X2 

fri = Cr2 = fro = C 



Ma 



(34) 



Since — fr/o is required, it indicates that there are 
some C = 2 band which can never support any PCI state 
with C4 symmetry, at least not one constructed by WSR. 



For example, consider the model: 

'H-{kx, ky) =(sin^ — sin^ ky)ax + 2sinA;a; sin kyay 

+ {1 + cosk,j; + cosky)(Tz (35) 

with Rc^ = (Tz- The lower energy band has ^Po = 
— 1, fr/o = Ij so that no C4 invariant gauge choice can 
be found. 

In summary, for a C = 2 Chern band on the square lat- 
tice, a scheme of mapping it into two C — 1 bands while 
preserving the C4 rotation symmetry does NOT always 
exist. A possible obstruction can be detected by calcu- 
lating the eigenvalue of the C4 operation at high sym- 
metry points in the BZ. The question remains: what are 
the conditions that guarantee a C4 symmetry preserving 
splitting scheme of a Chern band with a generic Chern 
number C? It is very likely that a generic Chern number 
will require a enlargement of unit cell that is incompat- 
ible with the rotation symmetry and, thus, will forbid a 
symmetry preserving splitting of the Chern band. This 
may be a reason to believe that the topological nematic 
state is generic. 



V. CONCLUSION AND DISCUSSIONS 

In conclusion, we propose a scheme to select the gauge 
of the Bloch states of C = 1 Chern bands so that the 
mapping from LL to the Chern bands respects the C4 
rotation symmetry of the square lattice. Consequently, 
FCI states obtained through this mapping are invari- 
ant under the C4 rotation, and thus, are better trial 
wavefunctions to be compared with the numerics. Also, 
through the C4 invariant coherent states, a C4 symmetric 
pseudo-potential formalism can be constructed. Our dis- 
cussion mainly focuses on the C — 1 Chern bands with 
C4 symmetry on square lattices. Similar construction 
can be done for those with C'e symmetry on triangular 
or hexagonal lattices. The discussion is generalized to 
higher Chern number case. A C = 2 band can be re- 
duced to two C = 1 bands by WSR. The C4 symmetry 
property of the C — 2 band can be determined by that 
of the two reduced C — 1 bands. C4 symmetry imposes 
some constraints on the direction of the Wannier states 
in WSR, and also some constraints on the rotation eigen- 
values at the high symmetry points. The conditions that 
guarantee a rotationally symmetric splitting scheme for 
a Chern band with an arbitrary Chern number C needs 
further investigation. 
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